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In this paper the properties of dark solitons in superfiuid Fermi gases with spin-imbalance are 
studied by means of a recently developed effective field theory [S. N. Klimin et al., Eur. Phys. J. 
B 88, 122 (2015)] suitable to describe the BEC-BCS crossover in ultracold gases in an extended 
range of temperatures as compared to the usual Ginzburg-Landau treatments. The spatial profiles 
for the total density and for the density of the excess-spin component, and the changes of their 
properties across the BEC-BCS crossover are examined in different conditions of temperature and 
imbalance. The presence of population imbalance is shown to strongly affect the structure of the 
soliton excitation by filling its core with unpaired atoms. This in turn influences the dynamical 
properties of the soliton since the additional particles in the core have to be dragged along thus 
altering the effective mass. 
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I. INTRODUCTION 

Solitons are nonlinear localized excitations that prop¬ 
agate trough a medium without changing their shape: 
these properties make them extremely valuable in the 
study of the interplay between coherence and interactions 
and hence they were examined in an extremely broad set 
of physical systems, ranging from macroscopic solitary 
waves in water canals, to pulses in ratical fibers [I| , to 
excitations in conducting polymers [3- 
In the field of ultracold atoms a great deal of attention 
has been focused on dark solitons: these appear as lo¬ 
calized density dips propagating on a stable background 
and are characterized by a phase jump across the den¬ 
sity minimum. While the first experimental studies of 
dark solitons in BEG condensates date back to the first 
decade of the 2000s their investigation in fermionic 
ultracold gases is more recent. In 2013 an experiment 
at MIT Q detected long-lived solitary waves that were 
initially identified as dark solitons; a later paper by the 
same research group Q explained the surprising long life 
and high effective mass of the observed excitations by 
correctly reinterpreting them as solitonic vortices: prod¬ 
ucts of the decay of dark solitons via the snake instabil- 
iW @,[3, already predicted and experimentally detected 
[m in BEG condensates. The main technique employed 
to create solitons in ultracold gases is phase imprinting 
il. The starting point is an ultracold gas in which the 
interparticle interaction can be tuned by using a Eesh- 
bach resonance. A laser beam is shone on one half of the 
atomic cloud, thus locally changing the phase of the order 
parameter. The soliton is then observed via absorption 
imaging after a time-of-flight expansion [^-Q . In the case 
of Fermi superfluids, during the time-of-flight expansion 
a rapid ramp towards the BEG side of the resonance is 
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also employed 0, @|- In 2014 [I^ a technique was pro¬ 
posed in order to create stable dark solitons in fermionic 
systems by applying the phase imprinting method to just 
one of the two spin-components of the gas. 

From the theoretical point of view, solitons in BEC-BCS 
crossover superfluids have been widely studied by solv¬ 
ing the Bogoliubov-de Gennes equations (l3l - [l^ or - lim¬ 
ited to the deep BEG side of the crossover - the Gross- 
Pitaevskii equation [13, El]- In this article we are going 
to perform the study of solitons in imbalanced fermionic 
superfluids in the context of the effective field theory 
(EFT) presented in [I^ and already employed for the 
examination of the population-balanced case [2^. This 
theory has the advantage of being applicable to the whole 
range of temperatures below the critical temperature 
(correctly retrieving the results of the Ginzburg-Landau 
treatment in the vicinity of Tc [2l| and of the effective 
theories developed specifically for T = 0 [23 - [2l| ) and of 
being much less demanding than the BdG approach from 
a computational point of view. In this paper we study 
dark solitons in superfluid Fermi gases with the addi¬ 
tion of spin-imbalance, i.e. the possibility of having un¬ 
equal populations of spin-up and spin-down particles in 
the ultracold gas. Imbalance in ultracold Fermi systems 
was first experimentally engineered in 2006 but, 

because of its effects on the pairing mechanism, it has 
been a focus of major theoretical attention starting from 
the seminal works on critical fields in superconductors 

B [^ . up to more recent papers on ultracold fermions 
. Dark solitons, like vortices, are characterized by a 
dip in the order parameter’s profile: we will show that, 
in presence of population imbalance, this dip provides a 
favorable place where to accommodate the excess com¬ 
ponent particles. A travelling soliton has to transport 
these excess particles: we compute how this affects the 
basic dynamical quantities of the system. 

The present article is organized as follows: in Section HIl 
we give a brief sketch of the theoretical model employed 
to examine our system Section|nT]is dedicated to 
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the solution of the equations of motion and the definition 
of some of the relevant quantities that will be analyzed in 
the following; in SectionH^we study the various aspects 
related to the filling of the soliton in presence of imbal¬ 
ance and consider how the range of applicability of the 
EFT is affected by temperature and imbalance; finally 
the conclusions of our investigation are given in Section 
0 

II. MODEL 


with the fermionic Matsubara frequencies ujn = (2n + 
l)7r//3. The analytic expression for fi{/3,xX) and the 
recursion relation that links it to the other /s(/3, x, C) are 
given in (l^ . 

As already mentioned, the two main characteristics of 
a soliton are a sharp phase gradient and a dip in the 
spatial profile of the order parameter, therefore we need a 
formalism that can properly describe both aspects. This 
can be obtained by separating these two contributions as 


The formalism employed in this treatment is the one 
presented in [^ . Assuming slow variations of the su¬ 
perfluid order parameter in both time and space, the 
effective action for the system in the natural units h = 1, 
2m = 1, Ep = vp = hkp/m = 2, is 


SiP) 



+ H 


( 1 ) 


Moreover we can highlight the change in the amplitude 
by defining the modulus of the order parameter as its 
bulk value (that can be obtained from the solution of the 
gap equation) multiplied with the amplitude modulation 
function a as |dt(r,t)| = oo| a(r, t). From this it be¬ 
comes immediately clear that, for a localized excitation 
such as a soliton, air —>■ oo,t) = Ooo = 1- The real-time 
Lagrangian for the system is given by 


where /3 is the inverse temperature and the Hamiltonian 
T-L is given by 
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The first term appearing in this formula is the thermo¬ 
dynamic potential Hg that reads 


Hg = - 


dfc 
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log [2 cosh {PEk) + 2 cosh (/3C)] ’ 
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Here the dispersion relation for a free 

fermion, Ek = \/^k + I'I'P is the Bogoliubov excitation 
energy and Og is the s—wave scattering length that de¬ 
termines the strength and sign of the contact interaction. 
Since we want to consider a system with spin-imbalance 
we introduce different chemical potentials and for 
fermions with spin up and down respectively. These are 
in turn linked to the quantities /r and C by the relations 
M = (Mt + Mt)/2 and ( = - Mi)/2. 

The other coefficients appearing in m and @ are defined 
as 

/ dfr k'^ 

(4) 

/ dk 

^ (/3, a, C) - /i (/5, Afc, C)] , (5) 

/ dfr k'^ 

^-2ilUW,EU), ( 0 ) 

where the functions fs{P,xX) introduced in the last set 
of expressions are the solutions of 


fsiP,x,C) = 


[^(Wn ~ iC)^ + ■ 


(7) 


C = -K{a)a^—-n, (9) 

with the real-time Hamiltonian 

n = Hg(a) - . 

( 10 ) 

The coefficient k of the term with the time derivative, 
the quantum pressure Pqp and the superfluid density psf 
are related to the coefficients C, D and E (eqs. (0]), (O, 
(©) by 


K(a) =D |4'oo|^ , 

(11) 

Ps/(a) 1^1^ > 

(12) 

, , C'-4|4'p£’ , ,2 

(13) 

Pqp{a) = -" . 
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III. DARK SOLITON SOLUTION 


Imposing the condition 

a{x,t)—> a{x — vst), 9{x,t) — >0{x — vst) (14) 

requiring that the soliton propagates with constant ve¬ 
locity Vs along the direction x, and demanding that the 
total change of phase across the soliton 

A0 = 6{x —>• oo) — 9{x —>• —oo) (15) 

remains finite, we can solve the Lagrange equations for 
for a(x) and 6{x) as in The spatial dependence 
of the phase is described by 


6{x) = vs f dx' 

J —OO 


K{a{x'))a{x'Y — K(aoo) 
Psf{a{x')) 


n 


(16) 
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From the Lagrange equation for a we instead obtain the 
following relation that gives the position (i.e. the dis¬ 
tance from the soliton center) for each value of the am¬ 
plitude 


a: = ± 


V2. 


dab 


Pqp{a') 


X{a') - vlY{a') ' 


(17) 


1. Shape of the soliton 

Figure [TJa) shows the fermion density profile (normal¬ 
ized to the bulk density) of the soliton as a function of 
the distance from its center, for different values of the 
imbalance C- The fermion density is calculated within 
the mean-field local density approximation (LDA), using 
the formula 


In this expression the functions X (a) and Y (a) are de¬ 
fined as 


A^(a) — fls(a) fls(aoo) ^ 
[/c(a)a2 - Koo]^ 

Y(a) = 


2psf(a) 


(18) 

(19) 


and the amplitude at the center of the soliton oq = 
a{x = 0) is easily determined as the solution of X{ao) — 
u|y(ao) = 0. 

At a later stage of the present paper we are going to 
briefly consider a few aspects related to the dynamics of 
the soliton, therefore we give here the definitions for the 
soliton’s momentum and energy. 


-n / \ o /o A a/ pqp{a)Y{a) 

rs{vs) = 2V2vs do ==^ 

Jan ^X(a) — VgY (a) 


- TTKn 


Ssivs) = 2 V 2 f da 

Jao \/X(a) — VgY {a) 


( 20 ) 

( 21 ) 


For a detailed derivation of the above expressions we ad¬ 
dress again the reader to [ 2 ^. 

Before presenting our results it is necessary to make a 
further consideration. Since one fundamental element in 
the derivation of the EFT that we are employing is a gra¬ 
dient expansion of the pairing field up to second order in 
spatial and time gradients, we have to carefully exam¬ 
ine where we must consider the coordinate dependence 
of the coefficients C, D and E and where doing so would 
lead us beyond the limits of our approximation. From (HD 
and ([2]) it follows that we need to keep the coordinate de¬ 
pendence in the coefficient D and in the thermodynamic 
potential fls- On the other hand, for what concerns the 
coefficients C and A, it is fair to keep them equal to their 
bulk values C'(aoo) and E(aoo)- 


IV. RESULTS 

In this section we present the numerical results for the 
soliton in the spin-polarized Fermi gas and study the ef¬ 
fect of population imbalance on its properties. Moreover 
the range of applicability of the present effective field the¬ 
ory is determined and its dependence on temperature, 
imbalance and interaction is analyzed. 


n(LDA)=_^, (22) 

op 

In ( 2 ^ the non-zero value of the density at the center of 
the soliton - where, for low values of vs, the order pa¬ 
rameter approaches zero - was explained by inferring that 
the soliton gets filled by unpaired fermions. Consistently 
with this hypothesis we observe that the presence of im¬ 
balance enhances the value of the density at the soliton 
center. More in general, for increasing values of (C, we can 
see that the soliton gets filled with a growing amount of 
particles, becoming less deep and slightly broader. This 
last effect is shown more clearly in the inset of Fig. [TJa), 
which depicts the behavior of the inverse soliton width 
(taken at half the height of the density dip) While 
the density profile has been plotted only for the BEC side 
of the crossover, the inset shows that similar effects take 
place at unitarity and in the BCS regime. The values 
of the soliton width in the BCS regime are noticeably 
higher than those in the other cases. 

Figures [Il(b)-[T](c) also show the fermion density profile, 
but for different values of velocity vs and temperature 
T/Tp respectively. The results for the temperature de¬ 
pendence are shown in the BCS regime where the effect 
of temperature is more evident. The above-mentioned 
filling of the soliton for an increasing spin imbalance can 
be characterized more distinctly by calculating the den¬ 
sity difference 5n{x) between the “spin-up” and “spin- 
down” populations along the soliton dip. In figure UKa) 
we display the normalized Sn(x)/6n(0) profile for differ¬ 
ent values of the imbalance C- The left inset clearly shows 
how the relative density difference at the center increases 
monotonously with the imbalance in all different regimes 
of the BEC-BCS crossover. 

As the imbalance in the Fermi gas increases, so does the 
amount of unpaired particles that cannot participate in 
the superfluid state of condensated pairs. At higher tem¬ 
peratures, some of these normal state particles coexist 
with the condensate as a thermal gas, but any remain¬ 
ing excess of the majority component has to be spatially 
removed from the pair condensate. The soliton dip then 
turns out to be a suitable location to accommodate nor¬ 
mal state particles, and consequently fills up with an in¬ 
creasing amount of unpaired particles as the imbalance 
gets higher. As it is shown in the right inset, the width 
of the Sn{x)/6n{0) curves increases with ( too, in agree¬ 
ment with the earlier observed broadening of the soliton 
density dip. 

Similar to the density profile plots, we have also included 
figures of 5n{x) for different values of velocity vs and 
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temperature T/Tp, presented in figure [U^b) and[2jc) re¬ 
spectively. In this regard it is interesting to notice that 
the decrease in temperature leads to an increasing degree 
of localization for the distribution of excess-spin compo¬ 
nent particles. 

From the insets of Fig[TJb) and FiglJl^b) concerning the 
width of the soliton and ^sn of the excess component 
distribution, it appears that a critical value of the soliton 
velocity can be determined above which the soliton can¬ 
not exist (the width of the profiles n{x) and 5n{x) goes 
to infinity). The effect of imbalance on this quantity has 
been analyzed and the results are shown in Fig [3] again 
all across the BEC-BCS crossover. From the comparison 
with the behavior of the mean field bulk value of the order 
parameter in the same regimes (inset) it becomes clear 
that the value of the imbalance for which the critical ve¬ 
locity goes to zero is the same for which the 

minimum of the free energy corresponds to I'I'ool = 0 i.e. 
when the normal state becomes energetically favorable 
over the superfluid one. 


2. Validity range of the treatment 

As mentioned above, the theoretical model employed 
in the present discussion is based on the assumption that 
the profile of the order parameter changes slowly in both 
time and space, corresponding to a long-wavelength ap¬ 
proximation. In order to verify whether or not this as¬ 
sumption holds for the system under consideration, we 
compare the characteristic size of the soliton with the 
Cooper-pair correlation length ^pair (often referred to as 
Pippard length [13). The comparison ultimately enables 
us to determine the domain of validity of the EFT and 
identify it as the region of the {/3, C, (fcj’Os)”^} -space for 
which the size of the soliton is much larger than ^pair- 
To describe the characteristic size of the soliton we will 
use (and compare) two different definitions: the width 
at half height of the soliton dip (already used in the 
previous section) and the healing length ^phase- In this 
context, ^phase for the soliton at rest (vs = 0) has been 
calculated by modeling the spatial profile of the order 
parameter with a trial form of the amplitude modulation 

a(a;) = tanh ( —— | (23) 

\V2ivarJ 

and then minimizing the free energy with respect to 
the variational parameter As discussed in [3l| a 

rescaling coefficient has to be introduced in order to con¬ 
nect the variational parameter ^yar to the healing length 
^ph ase- To this purpose, the convention of rescaling ^yar 
to the value of ^phase at T = 0 in the BCS limit was 
adopted, obtaining the relation ^ptase = 1.175 ^„ar- 
The pair coherence length ^pair has been calculated in 
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FIG. 1. (Color online) (a) Density profiles at T = O.ITf, 
vs = 0.4 on the BEC side of the resonance (kFas)~^ = 0.5 
for different values of the imbalance parameter (b) Den¬ 
sity profiles at T = O.ITf, C = 0.15 on the BEC side of the 
resonance (IcfUs)”^ = 0.5 for different values of the soliton 
velocity, (c) Density profiles at (" = 0.15, vs = 0.4 on the 
BCS side of the resonance (/cfUs)”^ = —0.5 for different val¬ 
ues of the temperature. The insets show the behavior of the 
inverse soliton width (^n)~^ hr the BCS (blue dotdashed line), 
unitarity (green dashed line) and BEC (red line) regimes as 
a function of vs and T respectively. The position x and 
widths are given in units of the imbalance parameter 
( is in units of Ep and the velocities are in units of uf/2. 
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FIG. 2. (Color online) (a) Excess component density at 
T = O.ITf, Vs = 0.4 on the BEG side of the resonance 
{kFas)~^ = 0.5 for different values of the imbalance, (b) Ex¬ 
cess component density at T = O.ITf, C = 0.15 on the BEG 
side of the resonance = 0.5 for different values of the 

soliton velocity, (c) Excess component density at (" = 0.15, 
vs = 0.4 on the BGS side of the resonance (/cfOs)”^ = —0.5 
for different values of the temperature. The left (right) in¬ 
sets show the behavior of the ratio 5n{0)/n{0) (inverse width 
(Ci5n)~^ of excess component distribution) in the BGS 
(blue dotdashed line), unitarity (green dashed line) and BEG 
(red line) regimes as a function of the imbalance, soliton veloc¬ 
ity and temperature respectively. The position x and widths 
^Sn are given in units of kp^, the imbalance parameter (( is in 
units of Ef and the velocities are in units of vfI‘2. 



FIG. 3. (Color online) Critical velocity of the soliton as a 
function of the imbalance at T = O.ITf across the BEC- 
BCS crossover. The inset shows the corresponding behavior 
of l^tool- The imbalance parameter ^ is given in units of Tf, 
the velocities in units of wf/2 and |'l/oo| in units of Ef- 


[ 3 ^ in terms of the pair correlation function 

m(^) = - (I) (24) 

by following the definition 


^pair — 


I f drr^g-fi(r) 

fdrgu(r) 


(25) 


At mean-field level, in the context of our effective field 
theory, this quantity can be obtained as 


fdkk^ (4kaf2(/3,Ek,Of 
fdkk^ 


(26) 


In Figl4^ and Figldjo the behaviors of the inverse pair 
coherence length soliton width and healing 

length function of the interaction strength are 

compared for two values {T/T^ = 0.1 and T/Tc = 0.95 
respectively) of the ratio between the temperature of the 
system and the critical temperature. As expected a very 
good agreement is found between the variationally deter¬ 
mined healing length and the soliton width Only in 
the vicinity of the unitarity regime we see a sizeable dif¬ 
ference between these two quantities. This can be under¬ 
stood since, while in both the BCS a and BEG [13, [11 
limits theoretical predictions state that the soliton am¬ 
plitude profile is described by a hyperbolic tangent as in 
(1^51) . at unitarity and in the intermediate regimes this is 
generally not true. 

For low values of the temperature (Fig|4^), i~air ^^nd 
are very close in the unitarity and BCS regimes but de¬ 
viate from each other when going towards the BEG side 
of the resonance, therefore limiting the validity domain 
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of the EFT to this region. On the other hand, from the 
data sets relative to T = 0.95Tc (FigHjD) it is apparent 
that the soliton width remains noticeably higher than the 
Pippard length across the entire depicted window of the 
BEC-BCS crossover, guaranteeing the reliability of the 
predictions of the EFT in the whole domain. To have a 
better understanding of the validity range of the EFT, in 
FiglS]the ratio between the pair coherence length and the 
healing length is plotted as a function of both tempera¬ 
ture (normalized to T^) and interaction strength. The 
validity domain of the EFT can be intuitively identified 
with the dark blue/purple region in the contour plot. 

For what concerns the dependence of ^pair on C, for the 
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FIG. 5. (Color online) Contour plot depicting the behavior 
of the ratio ^pair/^phase between the pair coherence length 
and the variationally determined healing length as a func¬ 
tion of the inverse scattering length Ua and the temperature. 
The values of the scattering length span the entire BEC-BCS 
crossover, while the temperature range goes from 0 to the 
critical temperature Tc. 


Figini for certain values of T a non-monotonic behavior of 
Cpair(C) is observed. This can be explained by examining 
the relation between the order parameter and the imbal¬ 
ance: in these configurations the modulus of the order 
parameter varies slowly in the entire range of values of C 
with the exception of the region in the immediate vicin¬ 
ity of the critical value where it suddenly drops to 

zero. 


3. Dynamical properties 

For a balanced Fermi gas, it has been demonstrated 
that, within the effective field theory, a soliton obeys 
the energy-momentum relation of classic Hamiltonian dy¬ 
namics 


FIG. 4. (Color online) (a) Comparison between inverse pair 
coherence length and inverse soliton width (in units of kp) 
across the BEC-BCS crossover for T/Tc = 0.1. The full 
black line represents the inverse healing length without im¬ 
balance. The dashed blue line describes the inverse soliton 
width. The dashed green line represents the inverse pair co¬ 
herence length. The dotted red line represents the inverse pair 
coherence length with imbalance C, = 0.5E_f. (b) Compari¬ 
son between inverse pair coherence length and inverse soliton 
width across the BEC-BCS crossover for T/T^ = 0.95: the 
dashing/color code is the same as for (a). 

situations considered in FiglHwe see that the presence of 
a population imbalance modifies just weakly the behav¬ 
ior of the pair coherence length. However, as shown in 


d£s 

oVs 


(27) 


i.e. the soliton behaves as if it were a classic parti¬ 
cle. Extending the formalism and calculations to a spin- 
polarized gas, we numerically confirm that this relation 
is still satisfied in the presence of imbalance. With this in 
mind, the effective mass Ms of the soliton can be defined 
in relation to the soliton momentum Vs (1^01) and to the 
soliton energy £s m as in [T^ : 


OVs Vs OVs 


(28) 


The soliton’s motion can then be treated as the motion of 
a classical particle with effective mass Ms moving with 
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FIG. 6. (Color online) Dependence of the pair coherence 
length on the imbalance parameter ^ and temperature T at 
{kFas)~^ = —1 on the BCS side of the resonance. The non¬ 
monotonic behavior of ^pair(C) is apparent. Lengths are in 
units of and ^ is in units of Ef- 


velocity vs- 

In Fig[7]the behavior of the effective mass (which is a neg¬ 
ative quantity) calculated at vs = 0 is examined across 
the BEC-BCS crossover for different values of the tem¬ 
perature. The effect of temperature becomes more evi¬ 
dent in the interval of values of the interaction parameter 
ranging roughly between 0 and 1. For low temperatures 
{T/Tp = 0.01) we notice the appearence of a peak cen¬ 
tered around (fc^Os)”^ = 0.5 while for higher tempera¬ 
tures a monotonic increase of the absolute value of Ms is 
observed as we move from the BCS side towards the BEC 
side of the Feshbach resonance. The effect of imbalance 
on Ms is shown in FigjS] where the absolute value \Ms\ 
is plotted as a function of {kFas)~^ for various values of 
As for temperature, also in this case the change in the 
behavior due to C is more evident in the vicinity of the 
unitarity regime. While in the BEC and BCS limits we 
observe that \Ms\ is larger for small values of the imbal¬ 
ance and decreases with increasing (( as could be intu¬ 
itively expected in consideration of the discussion about 
the filling of the soliton by imbalance, in the intermediate 
region we see that the behavior is much less straightfor¬ 
ward. In the immediate vicinity of (fc^Us)”^ = 0 in fact, 
the higher value of \Ms\ is not reached for a balanced 
system but for a highly imbalanced one (C = 0.4). 

A final remark has to be made for what concerns the in¬ 
fluence of the imbalance on the period of oscillation of the 
solitonic wave in the trapped condensate. This period Tg 
can be calculated according to the formula in [9|: 


Ts 



where Tp = 2 'kIujt is the period associated with the 
frequency of the trap and Ns is the (negative) amount 
of particles in the soliton dip (per unit area), defined as 
/_^[n(a;) —n(oo)] dx . The quantity mNs represents the 
physical mass of the soliton, which in general is different 
from the effective mass. Since both the physical and ef¬ 
fective mass of the dark soliton are negative, their ratio 
in (1^^ will be positive. It must be remarked that our 
present treatment, because of the assumption da, does 
not allow for oscillations of the soliton in the (y, z) plane. 
Therefore the following calculations can be interpreted 
as a limiting case for a strongly confined {quasi-ID) gas 
hence not allowing the snake instability to develop. 
Figure [9] shows the ratio of the soliton oscillation fre¬ 
quency to the trap frequency for the different regimes 
of the BEC-BCS crossover. In the BEC and unitarity 
regimes, the oscillation period increases monotonously 
for increasing values of C. In the BCS regime, the 
changes in the oscillation frequency are very small com¬ 
pared to the other cases. As opposed to what is stated in 
[Tsl . [ 3 ^ Is^l - and observed in Q in relation to solitonic 
vortices -, from FigjS] we see an overall increase in the 
ratio lvs/ojt as the system goes from the BEC towards 
the BCS regime. However, at low temperatures the de¬ 
pendence of (jJs/iVT on the interaction parameter is not 
monotonic: the ratio is found to first decrease from the 
deep BCS regime until it reaches a minimum for a value 
of (/cfOs)”^ around 0.5 and then to increase towards the 
theoretically predicted value of l/-\/2 as (/cfOs)”^ reaches 
high positive values. A sign of this discrepancy can be 
observed by considering the graphs relative to the soli¬ 
ton energy £s in where the energy calculated in the 
framework of the present EET is compared with the re¬ 
sults obtained from the solution of the BdG equations. 
From this comparison it emerges that in both the BCS 
and unitarity regimes the slope of £s{vs) according to our 
theory is less steep than the BdG one. As it is clear from 
(1^51) . this translates in an underestimation of the effective 
mass Ms and, in turn, in an overestimation of the ratio 
ujs/ujt in these regimes. To understand the observed dif¬ 
ference between the predictions of the present EFT and 
the other results found in literature, three more elements 
have to be taken into account. First it has to be stressed 
that equation (1^ is obtained by considering a soliton 
oscillating inside a trap, a condition that does not match 
our assumption p4l) requiring the excitation to move at 
constant speed vs along the x direction. Also, due to the 
bosonic nature of the EFT, Andreev bound states cannot 
be straightforwardly included in the treatment: however 
their presence is predicted to give a sizeable contribution 
to the effective mass of the soliton Finally in 

the present treatment we have fixed the value of the back¬ 
ground density far away from the soliton Uoo = I/Stt^: 
for a trapped system a correction should in principle be 
included in the calculation of Ns to account for the vari¬ 
ation of this quantity across the BEC-BCS crossover. 


(29) 








FIG. 7. (Color online) Absolute value of the effective mass 
of the soliton as a function of the interaction parameter 
for different values of the temperature without im¬ 
balance = 0). The mass is given in units of 2m. 



(^F^s) ^ 

FIG. 8. (Color online) Absolute value of the effective mass 
of the soliton as a function of the interaction parameter 
(IcfUs)”^ at temperature T = O.ITf for different values of 
the imbalance (given in units of Ef)- The mass is given in 
units of 2m. 


V. CONCLUSIONS 

In this paper we have analyzed the effect of tempera¬ 
ture and population imbalance on the properties of dark 
solitons in ultracold Fermi gases across the BEC-BCS 
crossover. The study was performed in the context of 
an effective field theory illustrated in based on the 
only assumption of a slow-varying pairing field. The fact 
that no hypothesis was made requiring the pair field 
to be small, in principle extends the range of validity 
of this theory with respect to the widely employed 
Ginzburg-Landau and BdG approaches enabling us to 
consider also the effect of temperature on the system. 



FIG. 9. (Color online) Ratio of the soliton oscillation fre¬ 
quency ujs to trapping frequency ujt, calculated as uis/^t = 
mNg/Ms as a function of the imbalance ^ (in units of Ef) 
at T — O.ITf and vs = 0.6 across the BEC-BCS crossover. 
Velocities are in units of 2vf. 


To test the intrinsic reliability of the predictions of the 
EFT on dark solitons, a comparison has been made 
between the characteristic size of the excitation and 
the pair coherence length [3l|, , resulting in a better 

understanding of how temperature and interaction affect 
the domain of applicability of the theory. In particular, 
while the calculations relative to the BEG regime are 
valid for all values of the temperature from 0 to Tc, at 
unitarity and on the BGS side of the resonance this is 
true just for temperatures close to the critical one. 
Based on analytic expressions for the amplitude and 
phase spatial profiles, the density and the density 
of the excess-spin component were obtained in LDA 
approximation using for the bulk value of the order 
parameter the mean field results. By systematically 
analyzing the density profiles we have observed how 
increasing the imbalance (and consequently decreasing 
the number of particles available for pairing) results 
in a filling of the soliton core that thus proves to be a 
convenient place where the unpaired particles can be 
stored. This translates into a decrease of the modulus 
of the effective mass of the excitation with increasing 
imbalance. However, in the crossover region in the 
vicinity of the unitarity regime we observe that the 
effect of the imbalance on Mg is reversed. A discrepancy 
was observed between our predictions concerning the 
change in the soliton-to-trap period across the BEG-BCS 
crossover and those reported in other papers all based 
on the solution of the time-dependent Bogoliubov-de 
Gennes equations [H, HI, [s^ at zero temperature. 
Keeping in mind the experimental setup employed in the 
investigation of solitons in ultracold gases a substantial 
part of our results was presented as a function of the 
soliton velocity vs so to have the possibility of a direct 
comparison with future experimental results. 
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